Stability and The Existence of Coherent Structure in Demixed State of Binary BEC 
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From a linear stability analysis of the Gross Pitaevskii equation for binary Bose Einstein conden- 
sates, it is found that the uniform state becomes unstable to a periodic perturbation of wave number 
k ii k exceeds a critical value kc- However we find that a stationary spatially periodic state does not 
exist. We show the existence of pulse type solutions, when the pulse structure for one condensate 
is strongly influenced by the presence of the other condensate. 

PACS numbers: 



After the realization of Bose Einstein condensation 
(BEC) in dilute atomic gases of Rb [H, |1|, Na Q, Li 
3 , there has been a lot of interest in the field of binary 
mixtures of the condensates. The physics of trapped di- 
lute condensates which are inhomogeneous, finite sized 
systems has become one of the most exciting fields of 
research in the last decade. The first two component 
condensate was produced with different hyperfine states 
of i?6®^. The experiments on i?6®^ atoms occupying the 
hyperfine states \F = l,mp — —1) and \F = 2, to/ = 1) 
have been done by JILA group(Myatt et al.,1997) [Ij. 
Later the group in MIT has been able to trap atoms 
in F=l with different mp values using optical methods. 
The repulsive interaction between the hyperfine states 
was reported and since then the stability of binary Bose- 
Einstein Condensates [3,Q has been extensively studied. 
Specially in case of two species BEC there is type of new 
phenomenology associated with it. Lots of theoretical 
and experimental works have been done on Na- Rb mix- 
ture [7| and on K-Rb mixtures [ll|, [I4I . Infact there is a 
range of interspecies interaction strength in which a Na- 
Rb mixture can be stable in harmonic trap [8|. In case 
of demixed state of binary BEC there is a tendency of 
spatial separation of the condensates, the interatomic in- 
teraction playing the vital role. The possibility of phase 
separation of the components and the densitypatterns of 
the system has been theoretically reported [3, [13] • The 
spatial separation between the condensates has been ver- 
ified experimentally by the JILA experiment [l3|. Re- 
cently the miscibility and immiscibility of the dual species 
{Rb^^ and Rb^'^) BEC depending on the strength of the 
external magnetic field and hence on the scattering length 
has been reported experimentally IJ] . In this paper we 
are able to show that under the condition A < (the 
same condition as experimentally observed [IJ) the phase 
separation occurs in dual species BEC and there is lo- 
calization of one species in presence of other in a pulse 



structure, where A = ^^ — 1. We also have depicted 

how the density of one species affects the localization of 
other species. 

In many physical systems the spreading of wavepack- 
ets occurs in the dispersive(linear) medium because dif- 
ferent component of wavepackets have different veloci- 
ties. It is the speciality of nonlinear systems that there 
exist some coherent structures [15;]. They are the spe- 
cial type of solutions of pattern forming system. Pulse, 
front, source and sink are the four basic things to deal 
with in the study of coherent structures. Several solu- 
tions like domain walls [la], antidark or grey solitons of 
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one component bound to a bright or dark solitons [17] 
in the other, have been found for two component Gross 
Pitaevskii equation. A particularly interesting variant 
involves introducing a linear coupling between the two 
condensates. In the case of this coupling being given a 
time dependence, it has been shown that one can induce 
a Rabi switch between the two condensates. The pulse 
type solutions are extremely important as they behave 
like the solitons and propagate within the medium keep- 
ing their shape invariant. So if there exists a pulse struc- 
ture, the multicomponent system can have a species wise 
localization. We have found both analytically and nu- 
merically pulse type solutions satisfying Gross Pitaevskii 
(GP) equation. 

When the atoms are cooled to extremely low temper- 
atures the de Broglie wave length associated with the 
atoms starts to increase. At temperatures below Tc (crit- 
ical temperature) the wave functions start to overlap. 
Under these conditions the atoms start to condense in a 
single quantum state which is known to be the Bose Ein- 
stein condensation. The nonlinear interaction between 
the atoms in this regime is due to the s-wave scattering 
between the atoms and the interactions are described 
by the single parameter a, called the s-wave scattering 
length as the interaction takes place at a very low ki- 
netic energy. The interaction plays a vital role. For a 
given atomic species under attractive interaction a col- 
lapse may occur if the particle number in the condensate 
exceeds a certain critical value Nc i.e, at high density 



state though for confined atomic gases the BEC is real- 
ized and also the collapsing dynamics has been observed 
[Gerton, et al. 2000]. Most experiments have been done 
with the atoms in repulsive interatomic interaction where 
the Thomas-Fermi approximation holds good. Here we 
restrict ourselves to a repulsive interaction. 

The GP equation for describing the binary mixture is 
written in terms of two fields ^o(r,i) and '4>i{r,t) which 
are complex valued functions representating the expec- 
tation values of the corresponding quantum fields. In 
absence of any external potential coupled GP equation 
has the following well known form, 



ihdt^Po = (-|^v2+go|V'o|'+5oi|^i|')^o (1) 



where, 



,91 = and goi = 

mm m 



go = 



ao and ai are the scattering length of the particles be- 
longing to state and state 1. aoi is determined when 
an atom in the state is scattered by the atom in the 
state 1. The validity of GP equation depends on the fact 
that s-wave scattering length is much smaller than the 
average interparticle separation. 

A spatially uniform solution has the form: 



'Mjt 1 
V'j = e " -7y ; J = 0, 1 

The chemical potentials satisfy, 

A*o =.9o|V'oP +3011-01^ = {9o+9ai)/V 

A^i == e + gi|^iP +5011^^0^ = {<^ + 9o+9oi)/V 



(3) 



(4) 



To test the stability of the uniform state, we introduce 
perturbation. 






J = 0, 1 



(5) 



Linearizing Eqs. (1) and (2) in dtpj and its complex 
conjugates, we get 
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Consistency requires, 
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(6) 



If oj^ turns out to be real then the mixture will be stable 
for all k. Hence we have to analyze whether w^ is real or 
not. It can be found from Eq. (6) that for all values of k 
the mixture is not stable. The mixture is stable only for 
those values of k for which 



2K 

K'^ + —{90 + .91) < (ffoi - .9051) 



(7) 



So, by taking constant (j) solution the stability of the mix- 
ture is not possible for all values of k. Thus we need to 
Analise the solution having spatially dependent terms, to 
see if there is a spatially periodic stationary state. 

Accordingly, we try the stationary spatially periodic 
solutions. 



ipo- 



1 



— = + 64>o^i{x) 



e "^ 



(8) 



Where 90o,i(^) are spatially periodic functions having 
the expansions. 



S(t)o{x) — y^ An cos{nkx) 
6(j)i{x) — y ^ Bn cos{nkx) 



It should be clear that the wavenumber k would be found 
from Eq. (6) with lo = and thus immediately leads to 
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This can never yield two positive values for k^ and hence 
does not support a stationary spatially periodic state. 
Consequently we seek a solution which is both space and 
time dependent and is uniformly propagating. 

We have taken the coupling constants to be -|-ve imply- 
ing a repulsive interaction between the atoms. But there 
may be a regime where the strength of repulsive inter- 
actions are too high to have the spatial inhomogeneity 
within the system. This results in the demixing instabil- 
ity within the system. Condensate and 1 can separate 
spatially. It is found previously that spatially uniform 
condensate is stable when goi ^ (5o<7i)^ • But when 
this condition is not satisfied there is a tendency of spa- 
tial separation of the condensates and the resulting state 
is called the demixed state. The tendency of separation 
of the condensates has been found in experiment by JILA 
group [l3J . 

Now we investigate the possible structure for such a 
demixed state. For this we have taken. 
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with 



with 



MO = A{Oe'^o(€ 



(11) 
(12) 



where (x — vt) — ^. i.e, we are considering two conden- 
sates moving with same velocity and in same direction. 
Inserting the above equations in Eq. (1) and (2), the flow 
equations for the binary system becomes, 
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d^qi = Q 
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(17) 
(18) 



where, d^(pj — qj for j=0,l and kq = ^d^A, ki — -^d^B. 
Inserting the value of kq in Eq. (15) we get the value of 



<?o = 



Where, 



and from Eq. (14) 
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After a little bit calculation from Eq. (19), the relation 
comes out to be 



= ra 
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Where, 
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Now we can write the pulse type solution for conden- 
sate which will be of the form of, 



A^ = (Ml + Ma coshA/gO"^ 



(25) 



Ml = -77^ 
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(27) 
(28) 



The structure of this pulse of one condensate thus de- 
pends on the amplitude of the other pulse at that location 
and at that time. As the pulses are like the solitons of the 
system so they will propagate keeping their shape con- 
stant for long time. Hence both of the condensates will 
localize and the localization will be gaussian in shape. 

A careful observation of Eq. (23) reveals the fact that 
there is a competition among the terms. Since the second 
term in the parenthesis of the expression of tq is the 
kinetic energy of the atoms, we can neglect it compared 
to others in the low temperature regime. Hence rg comes 
out to be {2goiB^ — 2^q). Now three cases may happen 
according to the competition between the terms. 

Case 1: When /io — goiB^, Mi and AI2 turn out to be 
00, i.e, according to the Eq. (25) A'^ is zero. So the pulse 
structure vanishes. 

Case 2: When ^0 < goiB^, Mi is -ve. Since M2 = 
±Mi, for Ml < and M2 < 0, A^ turns put to be -ve 
which can't happen and for Afi < and M2 > there is 
a infinite divergence at ^ = 0. 

Case 3: Only for the remaining condition /ig > goiB^ 
the pulse structure exist. Now we will express the con- 
dition in more well known form which is the evidence 
of the spatial separation of the two species being pulse 
structured. 

For the condensate with density uq and condensate 
1 with density ni 
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Now the order parameter tpQ — \/No(po and V'l — VNi(j)i . 
For the constant solution i.e., for = -=, |^oP = no = 

\A\^ = A^ (A being real) and |V'iP = ni = \B\^ = B^{B 
being real). Multiplying Eqs. (29) and (30) now the 
condition becomes ^qi < 5o5i which implies A < 0, 
i.e., the condition for which immiscibility between the 
species [14] because of strong interspecies interaction oc- 
curs and hence there exists the pulse type localization of 
each species. 

Scaling B' = 1 - Stt x 10^21^2 and B = 2 x lO^S' we 
have plotted the pulse structure (Eq. (25)) for different 
values of B which is the scaled version of B. In figures A 
and B is represented by A and B. 

Fig. 1 depicts the structure of the pulse of conden- 
sate changing with the amplitude of condensate 1 . As 




FIG. 1: Sensitivity of the pulse structure on the values of B. 



Taking, m=23 amu.. 



10"' m/s, /x ~ KbT, T=100 nK, 



ao,aoi ~ 10 m. We have plotted A with different values 



of B, ^ : 



Ott X lO'^A. 



the value of B increases the structure of the pulse be- 
comes more and more steep i.e, when the amplitude of 
condensate 1 is high it will make the atoms of the other 
condensate to localize in space while when the amplitude 
of condensate 1 is low the atoms of condensate has 
enough space to spread over. So the localization of one 
condensate is greatly affected by the presence of another 
condensate. 

In Fig. 2 it is obvious that at the lower value of B i.e, 
at B=0.2, 0.6 the atoms of condensate spreads over 
the space. When the value of B increases the spreading 
decreases and the atoms get gradually localized. At 
B=1.8, 2.2 the localization is very prominent. Now 
instead of condensate if we start with the GP equation 
for condensate 1 the pulse would generate for condensate 
1 and the qualitative nature of the pulse for condensate 
1 will remain the same. 

So far we are dealing with the condensates having 
the same velocity and moving in the same direction. 
Now we will investigate the situation when two conden- 
sates having different velocity propagates in opposite 
direction. We take condensate and 1 have the velocity 
vo and vi respectively and they are going in opposite 
direction. Here we have not encountered the case (con- 
densates approaching each other) of collision between 
the condensates. We will take the case when from a 
mixture of two condensates they are starting to go apart. 
So to deal with the case, we will take 



ipo 



'1^0' 



" Mx- vot) and ?Ai = e" 



""' fi{x + vit) 



Where, 



mo) 
/i(6) 






(31) 

(32) 




FIG. 2: 2D visualization of the pulse structure for different 
values of B. From upper left corner to lower right corner, the 
values of B are 0.2, 0.6, 1.0, 1.4, 1.8 and 2.2 respectively. 



After doing little bit of calculation, we will find 

r-olo = ^[-2^*0-^1-2 + 2.901 [5(^1)]'] (33) 
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In the expression for tq for both the condensates 2nd 
term is always -l-ve, where as 1st and 3rd term is always 
-ve. Hence for a particular value of vq and vi there will 
be the situation when tq will become zero and the pulse 
for the corresponding condensates will cease to exist. We 



have found that at 



-(goi^ — /^o) the pulse for 



condensate will not exist but that of condensate 1 re- 
mains and at v\ = —{gmA^ — Mi) pulse of condensate 1 
vanishes though that of condensate survives, vq and vi 
being of the order of 10~^ meter/sec. 

We have studied the stability of binary BEG and we 
have found that for spatial independent solution, the 
system is not stable for all values of k. The system 
is stable only for those k which satisfy the condition 
K'^ + ^{90+ 9i) < (5oi - 9091) ■ By taking uniformly 
propagating solution of GP equation. We have observed 



the existence of pulse in the demixed state when the con- 
dition of immiscibihty holds. Under such condition the 
density of one species affects the localization of others. 
When two condensates moving in the opposite direction, 
at the low velocity both of the pulses survive and at high 
velocity regime, at a particular velocity (either vq for con- 



densate and vi for condensate l)one of the pulses ceases 
to exist. 

It is a pleasure to thank A. Bhowmick for helping in 
numerical section. Financial support by the S. N. Bose 
National Centre for Basic Sciences is gratefully acknowl- 
edged. 



[1] M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. 

Wieman and E.A. Cornell, Science, 269 198 1995. 
[2] C. J. Myatt, E. A. Burt, R. W. Christ, E. A. Cornell, [13] 

C.E. Wieman, Phys. Rev. Lett. 78, 586 (1997). 
[3] K. B. Davis, M. O. Mewes, M. R. Andrews, N. J. van [14] 

Druten, D. S. Durfee, D. M. Kurn and W. Ketterle, 

Phys. Rev. Lett. 75, 3969 (1995) . [15] 

[4] C. C. Bradley, C. A. Sackett, J. J. ToUett and R. G. 

Hulet, Phys. Rev. Lett. 75, 1687 (1995). [16] 

[5] R. Kaiser, C. Westbrook, and F. David, Coherent atomic 

matter waves, NATO Advanced Study Institute. [17] 

[6] P. Ao and S. T. Chui, J. Phys. B: At. Mol. Opt. Phys. 

33 535-544 2000. 
[7] R. Ejnisman, H. Pu, Y. E. Young and N. P. Bigclow, [18] 

Optical Express 330 2 1998 8. 
[8] C. K. Law, H. Pu, N. P. Bigelow, and J. H. Ebcrly, Phys. [19] 

Rev. Lett. 79, 3105 (1997) 
[9] F. Riboh and M. Modugno, Phys. Rev. A 65, 063614 [20] 

(2002) 
[10] D. M. Jezek and P. Capuzzi, Phys. Rev. A 66, 015602 [21] 

(2002) 
[11] G. Modugno, M. Modugno, F. Riboli, G. Roati, and M. [22] 

Inguscio, Phys. Rev. Lett. 89, 190404 (2002) 
[12] G. Thalhammer, C. Barontini, L. De Sarlo, J. Catani, F. 



Minardi, and M. Inguscio Phys. Rev. Lett. 100, 210402 
(2008) 

D.S. Hall, M.R. Matthews, J.R. Ensher, C.E. Wieman 
and E.A. Cornell, Phys. Rev. Lett. 81, 1539 (1998). 
S. B. Papp, J. M. Pino and C. E. Wieman, Phys. Rev. 
Lett. 101, 040402 (2008). 

W. V. Saarloos and P. C. Hohenberg, Physica D 56, 303- 
367 (1992). 

B. Deconinck, P. G. Kevrekidis, H. E. Nistazakis, and D. 
J. Frantzeskakis, Phys. Rev. A 70, 063605 (2004). 
P. G. Kevrekidis, H.E. nistazakis, D.J. Frantzeskakis, B. 
A. Malomed and R. Carretero-Conzalez, Eur. Phys. J. D 
28 181-185 2004 

O. Descalzi, J. Cisternas, P.Gutierrez, and H.R. Brand, 
Eur. Phys. J. S, p (e)cial Topics 146 63-70 2007. 
M. V. Hecke, C. Storm, W. V. Saarloos, Physica D 134, 
1-47 (1999). 

M. Inguscio, S. Stringari and C. E. Wieman, Interna- 
tional School of Physics 

F. Dalvano, S. Giorgini, Lev. P. Pitaevskii and S. 
Stringari, Rev. Mod. Phys. 71, 3 (1999) 
E. Timmermans, Phys. Rev. Lett. 81, 26 (1998) 



